Association of ultracold double-species bosonic molecules 
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We report on the creation of heterospecies bosonic molecules, associated from an ultracold Bose- 
Bose mixture of and *^Rb, by using a resonantly modulated magnetic field close to two Feshbach 
resonances. We measure the binding energy of the weakly bound molecular states versus the Fes- 
hbach field and compare our results to theoretical predictions. We observe the broadening and 
asymmetry of the association spectrum due to thermal distribution of the atoms, and a frequency 
shift occurring when the binding energy depends nonlinearly on the Feshbach field. A simple model 
is developed to quantitatively describe the association process. Our work marks an important step 
forward in the experimental route towards Bose-Einstein condensates of dipolar molecules. 

PACS numbers: 34.20.Cf, 34.50.-s, 67.85.-d 



Ultracold polar molecules hold the promise of a revolu- 
tion in the domain of quantum degenerate gases and pre- 
cision measurements. Degenerate molecules are sought 
primarily to produce a gas with strong long-range interac- 
tions, stemming from the coupling of electric dipole mo- 
ments that heterospecies dimers feature. Such molecules 
would create strongly correlated systems with a wealth of 
quantum phases fl| , provide candidate qubits [3, jsl , allow 
for a new generation of dipolar Bose-Einstein condensates 
(BECs) (jL and help in the search of the electron dipole 
moment [Sj. Starting from ultracold atoms, molecules 
have been successfully created following two different ap- 
proaches: photoassociation Q and magnetoassociation 
[7| , but few experiments have hitherto reported the pro- 
duction of Feshbach heteronuclear dimers. Two groups 
have reported the creation of fermionic KRb molecules 
[1, 0] , while the only bosonic dimer so far associated is 
ssRb^^Rb 0, which can not be dipolar since the con- 
stituents share the same electronic configuration. Het- 
erospecies bosonic dimers, i.e., the constituents of the 
dipolar BEG envisioned in Ref. [4|, have instead eluded 
experimental realization so far. 

A Bose-Bose mixture is particularly suitable to asso- 
ciate such dimers due to the high phase-space densities 
achievable, while the atom-dimer relaxation, which lim- 
its the lifetime of the molecules, can be strongly sup- 
pressed in optical lattices with a single atom pair per 
lattice site 11 ll. R ecent progresses in molecular stabiliza- 
tion schemes [12'| make the Bose-Bose ''^K^^Rb mixture 
truly promising for the experimental observation of BEGs 
of dipolar molecules. Following a different route, other 
experiments have very recently obtained ultracold het- 
erospecies Fermi-Fermi mixtures 13|, |14| that also can 



provide a way to compound bosonic dimers. 

In this Letter, we report on the production of het- 
erospecies "^^K^^Rb bosonic molecules starting from an 
ultracold mixture. In proximity of Feshbach resonances 
(FR's) at moderate magnetic fields, by adding a modu- 



lation to the Feshbach field [15| , we have converted up to 
12 000 '^^K^'^Wo pairs into dimers, i.e., 40% of the minor- 
ity ^^K atoms, at temperatures between 200 and 600 nK. 
We estimate the molecular lifetime to be at least 60 /is at 
atomic densities of ~ 5 x 10"'^^ cm~^ for each species. 

As described earlier [l^ > we start from a double-species 
MOT of ^^K and s'^Rb loaded by two separate 2D MOT's 
17| . A quadrupole magnetostatic trap moves the atoms 
to the center of a millimetric loffc trap, that confines 
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FIG. 1: (Color online): (a) Molecular radio-frequency asso- 
ciation: two free atoms with relative kinetic energy -Ekin are 
resonantly transferred to a molecular level by modulating the 
magnetic field. The inset shows the remaining atom num- 
ber after a variable association time. The solid line is an 
non-exponential fit to data (see text), (b) Spectrum for an 
association time of 30 ms, _B=78.25 G. The solid line is a fit 
with the lineshape described in text. 
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FIG. 2: (Color online): (a) KRb binding energy and 3-body 
losses for the Feshbach resonance at 79 G: triangles show the 
experimental binding energies, circles the atom number. The 
dash-dotted black line shows the theoretical prediction for 
i?b, the solid blue line is the same curve but translated by 
—0.24 G. The solid red line represents a Gaussian fit to the 
atom number. The vertical dashed lines show the position of 
the FR as determined by the 3-body losses (red, left) and the 
binding energy (blue, right), (b) Same as above, for the FR 
at 39 G. Note the widely different scale of a;-axis. Here the 
theoretical model for Eh is shifted by — 1 G. 



the sample during microwave evaporation on the *^Rb 
hyperfine transition and allows for sympathetic coohng 
between Rb and K. At a temperature of 1.5 /iK, we load 
an optical trap created by two orthogonal laser beams 
(A=1064nm, waist=110/im) and transfer *^Rb and ^^K 
from the |F = 2, m = 2) to the |1, 1) state by microwave 
and radio-frequency adiabatic passage, respectively, in 
presence of a 7 G bias magnetic field. All *^Rb atoms left 
in |2, 2) (^ 10%) are removed by a pulse of resonant light. 
No blast light is needed for "^^K atoms. At this stage we 
turn on a uniform magnetic Feshbach field of approxi- 
mately 78 G, close to the previously observed FR [ist . 
corresponding to a value of the interspecies scattering 
length a ~ 200 oq. By lowering the dipole trap power, we 
further evaporate the gas down to a temperature as low 
as 200 nK. At this temperature we have approximately 
4(3) X 10^ atoms of s'^Rb (^iR). 

To associate KRb molecules we bring the Feshbach 
field to the desired value below the resonance position Bq 
and modulate it by adding a radio-frequency (RF) field 
of typical amplitude ~130mG, as depicted in Fig. [IJa). 
This modulation pulse has a simple square envelope, with 



a typical duration of 20 ms. At magnetic fields B lower 
than the resonance value Bq , dimers are created when the 
modulation frequency /m is close to the binding energy 
E\^/hoi the molecular level. The production of molecules 
is revealed by the reduction of the total number of atoms 
N = A^K + A^Rb at the end of the RF pulse, as we scan the 
modulation frequency. A typical lineshape representing 
the molecular association peak is reported in Fig. [Tfb). 
Our molecules decay, most likely because of vibrational 
quenching induced by collisions with unpaired atoms Q , 
and leave the trap releasing the acquired binding energy. 

For several values of we measure the corresponding 
binding energy by fitting the experimental RF spectra 
with the lineshape model described later: the results are 
shown in Fig. [2] together with the theoretical predictions 
of the coUisional model [l3| ■ For both FR's the agreement 
is fairly good, provided the theoretical curve is translated 
by -0.24 G [Fig.IHa)] and -IG [Fig. [2tb)] along the B 
field axis. A precise determination of the position of the 
FR's is given by the points where the molecular energy 
curves fitting our data cross the i?b=0 threshold: 38.2 G 
and 78.67 G. To stress the sensitivity of this method, it 
is worth to remark that a deviation of 50 mG in the B 
field leads to a shift of ^50 kHz in the binding energy 
of the leftmost experimental point of Fig. [21^ a). In our 
experiment, the Feshbach field stability represents the 
main limitation on the precision of this method. We cali- 
brate the Feshbach field from the frequency of Rb hyper- 
fine transitions. Given the fluctuations observed during 
a single measurement and day-to-day, we associate an 
uncertainty of ±30 mG to the B field values. 

We compare these FR positions with the determina- 
tion obtained by the 3-body losses. We have repeated the 
loss measurements of Ref. |18| , the data are also shown in 
Fig. m For both FR's, the position obtained by the peak 
of the 3-body losses lays slightly below the B value ob- 
tained by the binding energy measurements. Actually, a 
more comprehensive analysis of the dynamics associated 
with 3-body losses is required to extract the FR posi- 
tion, especially for the the low-field broad FR. We also 
remind that the extrapolation to zero binding energy is 
performed under the assumption that the theoretical pre- 
dictions of Ref. [l9| are accurate except for a small shift 
in the B field. 

Our experimental RF spectra display several interest- 
ing and non-trivial features: (i) stark asymmetry and 
a pronounced broadening related to the finite tempera- 
ture of the atoms, (ii) a, shift of the resonant modulation 
frequency that increases with the modulation amplitude, 
and (in) additional association peaks at fractional fre- 
quencies of the binding energy. A precise description of 
the above mentioned features is crucial to accurately de- 
termine the binding energies from the measured RF spec- 
tra. 

The physical reason of the asymmetric lineshape shown 
in Fig. [TJ^b) is easily understood: as discussed later. 



3 



the resonant modulation frequency of a given atom 
pair depends on the pair kinetic energy i^kin, there- 
fore the spectrum is inhomogeneously broadened by the 
asymmetric Boltzmann distribution of kinetic energies 
VEkin exp(— i?kin/feT) Indeed for our typical tem- 
peratures, that is 200-600 nK, thermal broadening is the 
dominating contribution to the observed linewidth. The 
asymmetry went unobserved in related experiments of 
molecular association 0, [3] that employ degenerate or 
nearly degenerate gases. 

In the following we describe the simple model devel- 
oped to analyze our data by combining ideas from the 
theoretical description of photoassociation |213 with a 
model of the dynamics of molecule association in a har- 
monic oscillator well Our model is based on several 
simplifying assumptions. In the center-of-mass frame, 
we restrict to three quantum levels: one for each atomic 
species in the continuum and another one for the molec- 
ular bound state. Similar to [20| the multitude of contin- 
uum states is taken into account by the thermal average 
of the results. We start with a set of nonlinear differen- 
tial equations for the amplitudes aj {j = 1, 2) and m of 
the two atomic and the molecular state, respectively: 

icij — O cos(27r/int)a^m k^j (1) 
im — n cos{2Tr fnit)aia2 — {E{t)/h + i'j/2)m 

Here E{t) ~ Skin + E\^{t), ^'kin is the relative kinetic en- 
ergy of the atoms in their center-of-mass frame, Eh{t) = 
r][/^B + sin(27r/m<)]^ the time-dependent binding en- 
ergy at an average detuning AB from the FR. We use the 
quadratic form of E^, with curvature 77, valid in the uni- 
versal regime, and include the time variation due to the 
modulation of the magnetic field. For reasons of simplic- 
ity, we assume that the coupling strength 17, proportional 
to the modulation amplitude Bm, is the same for all states 
under consideration and that the molecule decay rate 7 
is independent of atomic and molecular densities. 

Numerical solutions of Eq. [1] describe several key fea- 
tures of our experimental observations: For sufficiently 
small modulation amplitudes this model predicts atom 
loss when the modulation frequency equals /o = (-Ekin + 
Eb)/h. The total atom number N{t) = \ai{t)\'^ + \a2{t)\'^ 
decays non-exponcntially: for balanced atomic popula- 
tions, we find 

in the limit <C 7. In our measurements we observe such 
an non-exponential decay, as shown in the inset of Fig.[T] 
This corresponds with the solution of the rate equation 
N = —K{fm)N'^ which is characteristic for molecule for- 
mation as a two-body process. The calculated atom loss 
coefficient K(/m) shows a Lorentzian dependence on the 
modulation frequency /m, with center and width given 
by /o and 7. In order to account for the thermal dis- 
tribution of the kinetic energies, we still fit our spectra 
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FIG. 3: (Color online): Lineshapes of the molecular associa- 
tion process for different values of temperature. Full circles 
represent data, taken after 30 ms of RF pulse length around 
5=78.2 G. Solid lines represents the fitted lineshape with the 
model described in the text: we take the binding energy Bb, 
the decay rate 7 and the couphng strength as free fit pa- 
rameters. The temperature is measured from the ballistic 
expansion of the clouds. Dashed lines show the fits assuming 
7 = 0. 

with Eq. [21 after replacing K(/m) by its convolution with 
the Boltzmann distribution: the results well describe our 
data, as shown by the fits in Fig. [3l 

Together with the binding energy, from fits we obtain 
an estimate for the molecular lifetime t = I/7 = 60 /is 
with an uncertainty of a factor 3. For comparison we 
report in Fig. [3] also the fit assuming no decay: it is 
evident that the finite 7 smoothens the left edge of the 
lineshape. Actually, from the spectra we can only infer a 
lower limit for the lifetime, that is r > 20 /us, but not an 
upper limit. Indeed the observed smoothening could be 
also explained by technical reasons, for example magnetic 
field fluctuations during the RF pulse. Instead, an upper 
limit T < 5 ms is deduced from the lack of an observable 
atom decay after the RF pulse is terminated. 

As mentioned above, other non-trivial features appear 
at high modulation amplitudes, that are also described 
by our model. First, we observe a shift of the modula- 
tion resonant frequency, increasing with the modulation 
amplitude (see Fig. [4]) , that can be understood in a very 
direct way considering the quadratic dependence of the 
binding energy on the magnetic field: the time-averaged 
value of the transition energy {E{t)) in Eq. [T] deviates 
from the value i?kin + ?/Ai?^ by an amount riB^/2. We 
directly verified this result by measuring the resonant 
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FIG. 4: (Color online): Shift of the resonant modulation fre- 
quency as a function of the modulation amplitude, measured 
at 5=78.30 G with 15 ms of RF pulse length and a tempera- 
ture of 350 uK. The solid line represents the parabolic fit. 



modulation frequency as a function of the modulation 
amplitude Bm (see Fig. As a consequence, for the 
high-field FR we extrapolate to zero modulation ampli- 
tude and correct the measured resonant frequencies by 
—6 kHz, corresponding to our modulation amplitude of 
0.13 G. For the low-field FR, no shift due to the ampli- 
tude modulation occurs since the molecular level is ap- 
proximately linear with the magnetic field. 

Moreover, we observe additional loss peaks occurring 
when the modulation frequency is at fractional values of 
the average transition energy /i/m = {E(t))/n with inte- 
ger n. This allows us to access a range of binding energies 
which lies outside the useable bandwidth of our excitation 
coil. The data points at large binding energy of Fig. [2] 
are deduced from such measurements with a modulation 
frequency corresponding to half the transition energy. 

Finally, we comment on the markedly different asso- 
ciation efficiency observed between the two FR's. On 
the low-field Feshbach resonance the dimers association 
shows a much poorer efficiency and requires an RF pulse 
of 1 s at maximum modulation amplitude 0.5 G. This can 
be understood considering that the coupling strength f2 
of Eq. [T] is significant only at field detunings AB where 
the bare states of the open and closed channel mix and 
the energy of the molecular level deviates from the linear 
dependence on the B field 23|| . As shown in Fig. [21 this is 
clearly the case only for the Feshbach resonance around 
79 G. Indeed we calculate that, in the range of our mea- 
surements, f2 is approximately a factor 10 smaller on the 
low-field FR, also taking into account the larger modu- 
lation amplitude employed. The reduction of association 
rate by two order of magnitudes is reflected in the need 
of an elongated RF pulse. 

In conclusion, we have reported the creation of het- 
erospecies ^^K^^Rb ultracold molecules, representing the 



first feasible starting point toward a BEG of dipolar 
molecules, since the different electronic structure of the 
constituents allows for a permanent electric dipole mo- 
ment in the vibrational ground state. On the most fa- 
vorable FR, with thermal clouds at 200 nK we converted 
up to 40% of the minority atoms. We carried out a de- 
tailed yet simple analysis of the association, highlighting 
a remarkable shift of the atoms-molecule transition that 
must be taken into account for a proper determination of 
the binding energy, when the latter depends nonlinearly 
on the Feshbach field. 

We plan to confine the atoms in an optical lattice 24 1 
to suppress the molecular decay due to inelastic collisions. 
This should grant enough time to drive the molecules to 
a lower ~ possibly the ground - vibrational states. Last 
but not least, molecular association will also serve as an 
important tool in the exploration of the quantum phases 
of the Bose-Bose mixture in the optical lattice. 

This work was supported by Ente CdR in Firenze, 
CNR under project EuroQUAM, EU under STREP 
NAME-QUAM, and INFN through project SQUAT- 
Super. We thank A. Simoni for sharing his theoretical 
results, and G. Varoquaux who contributed to the ex- 
periment, and our coUegues of the Quantum Degenerate 
Gases group at LENS for fruitful discussions. 

Note added. After the completion of this work, we have 
become aware of the production of bosonic molecules 
from the Fermi- Fermi mixture ^Li^'^K ^25i] . 
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